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CLr ABSTRACT 

2 ■ Context. The multi-wavelength study of old (>100 Myr) radio pulsars holds the key to understanding the long-term evolution 
of neutron stars, including the advanced stages of neutron star cooling and the evolution of the magnetosphere. Optical/UV 
^3 observations are particularly useful for such studies because they allow one to explore both thermal and non-thermal emission 
processes. In particular, studying the optical/UV emission constrains temperature of the bulk of the neutron star surface, too 
cold to be measured in X-ray observations. 

Aims. Aim of this work is to identify the optical counterpart of the very old (166 Myr) radio pulsar J0108— 1431. 
Methods. We have re-analyzed our original Very Large Telescope (VLT) observations (Mignani et al. 2003), where a very faint 
qq ' object was tentatively detected close to the radio position, near the edge of a field galaxy. 
I/-} , Results. We found that the backward extrapolation of the PSR J0108— 1431 proper motion recently measured by Chandra 
£\\ • (Pavlov et al. 2008) nicely fits the position of this object. Based on that, we propose it as a viable candidate for the optical 
counterpart to PSR J0108-1431. The object fluxes (U = 26.4±0.3; B » 27.9; V > 27.8) are consistent with a thermal spectrum 
with a brightness temperature of ~ 9 x 10 4 K (for R — 13 km at a distance of 130 pc), emitted from the bulk of the neutron 
star surface. 

Conclusions. New optical observations are required to confirm the optical identification of PSR J0108— 1431 and measure its 
spectrum. 
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1. Introduction to date. In particular, its position in the P-P diagram 

falls quite close to the pulsar death line. Indeed, with a 

The radio pulsar J0108-1431 was discovered by Tauris 40Q MHz luminosity of 0.15d 2 30 m J y kpc 2 , where d 130 is 

et al. (1994) during the Parkes Southern Pulsar Survey the distance in units of 130 pCj PSR J0108-1431 is the 

(Manchester et al. 1996). PSR J0108-1431 is certainly second faintest radio pulsar . 
one of the closest pulsars. Its dispersion measure of 2.38 pc 

cm -3 (D'Amico et al. 1998), the smallest known for a radio Soon after its discovery, there were several attempts 
pulsar, corresponds to a distance of 130 pc, according to to detect PSR J0108-1431 at other wavelengths. In X- 
the Galactic electron density model by Taylor & Cordes rays, the pulsar was not detected in the ROSAT All 
(1993) (180 pc in the model by Cordes & Lazio 2002). Sky Survey, while pointed UV, IR, and optical observa- 
The pulsar's period, P = 0.808 s, and period derivative, tions with EUVE (Korpela & Bowyer 1998), ISO (Koch- 
P = 7.44 x 10" 17 s s- 1 (D'Amico et al. 1998), correspond Miramond et al. 2002), and the SAO RAS 6m telescope 
to the spin-down age r = P/2P = 166 Myr, rotational (Kurt et al. 2000) did not reveal candidate counterparts, 
energy loss rate E = 5.8 x 10 30 ergs s~\ and surface Mignani et al. (2003) measured the most accurate radio- 
magnetic field B = 2.5 x 10 11 G. PSR J0108-1431 is interferometric position with ATCA and performed deep 
thus one of the oldest non-recycled radio pulsars known optical observations with the Very Large Telescope ( VLT) , 

but they did not claim a candidate counterpart. Recently, 

Send offprint requests to: R. P. Mignani; rm2@mssl.ucl.ac.uk PSR J0108-1431 was identified in a 30 ks exposure with 

* Based on observations collected at ESO, Paranal, under the Chandra ACIS detector (Pavlov et al. 2008), based on 

Programme 65.H-0400(A) its position close to the ATCA radio coordinates and on 



2 



Mignani et al.: Optical counterpart of PSR J0108-1431 




the high (> 300) X-ray-to-optical flux ratio deduced with 
the aid of the VLT images of Mignani et al. (2003). The 
pulsar's spectrum can be fitted with either a power-law 
with photon index T w 2.2 (for the hydrogen column den- 
sity iVn = 7.3 x 10 19 cm~ 2 estimated from the pulsar's dis- 
persion measure) or a blackbody with temperature kT 
0.28 keV and emitting area of ~ 50c?f 30 m 2 . The estimated 
0.3-8 keV luminosity, ~ 10 28 d 2 30 ergs s _1 , translates into 
the X-ray emitting efficiency rjx = Lx/E ~ 4 x 10~ 3 d 2 30 , 
somewhat higher than for most of other pulsars. The com- 
parison between the Chandra and ATCA coordinates of 
PSR J0108— 1431 yielded its first proper motion measure- 
ment (199±65 mas yr™ 1 ), corresponding to the transverse 
velocity of 129 ± 42 d 130 km s" 1 (Pavlov et al. 2008). 

In this paper we present a re-analysis of the VLT obser- 
vations of Mignani et al. (2003). The analysis and results 
are described in §2, while the discussion is presented in 
§3. 

2. Observations and data analysis 

Optical observations of the PSR J0108— 1431 field were 
performed in July and August 2000 with the VLT, using 
the original version of the FOcal Reducer/low dispersion 
Spectrograph (FORS1). At the time of the observations 
the instrument was equipped with a four-port 2048x2084 
CCD with a 0'.'2 pixel size and a field of view of 6'8 x 6'8 
in its standard resolution mode. Observations were per- 
formed in the Bessel U, B and V filters (see Mignani et 
al. 2003 for a detailed description of the observations and 
data reduction). As discussed by Mignani et al. (2003), 
the identification of the PSR J0108— 1431 counterpart was 
complicated by the fact that the ATCA radio position 
was close to the edge of a relatively bright elliptical field 
galaxy. A possible counterpart was noticed in the U and 
B bands only, « 0'.'3 northwest from to the radio position. 
However, the marginal detection significance of this object 
and its proximity to the galaxy made it difficult to estab- 
lish whether it was a genuine detection or a background 
enhancement, perhaps associated with the galaxy. For this 
reason, no identification with the pulsar was claimed. 

2.1. Astrometry 

The recent pulsar proper motion measurement by Pavlov 
et al. (2008) prompted us to re-investigate the possibil- 
ity that PSR J0108-1431 can be identified with the ob- 
ject marginally detected in the VLT data. We recomputed 
the astrometric calibration of our VLT/ FORS1 images us- 
ing as a reference the updated release of the GSC2 (ver. 
3.2), which provided a larger number of reference stars 
in the field and an improved coordinate accuracy, with a 
mean random radial error (for stellar objects) ctqsc — 0'.'3 
(Lasker et al. 2008). The detector coordinates of 20 se- 
lected GSC2 stars (all non-saturated, with a well-defined 
PSF, and evenly distributed in the field of view) were mea- 
sured by fitting a Gaussian function to their intensity pro- 
files using the Graphical Astronomy and Image Analysis 



Fig. 1. VLT U band image of the PSR J0108-1431 field 
obtained by Mignani et al. (2003). The pulsar candidate 
counterpart (U = 26.4 ± 0.3) is marked by an arrow. The 
cross shows the radio-interferometric position of Mignani 
et al. (2003) while the circle corresponds to the most re- 
cent pulsar position obtained with Chandra by Pavlov et 
al. (2008). The position of the candidate pulsar optical 
counterpart (epoch 2000.6) nicely fits the backward proper 
motion extrapolation. 

(GAIA) toofl The coordinate transformation between the 
detector and the celestial reference frame was then com- 
puted using the Starlink package ASTR0I"@. The rms of our 
astrometric fit was 0'.'14 for each of the two coordinates, 
corresponding to the radial error am = 0'.'21. The new as- 
trometric solution is consistent with the one obtained by 
Mignani et al. (2003), but it is a factor of 2 more precise, 
thanks to the larger number of GSC2 stars and to the im- 
proved astrometric accuracy. We estimated the overall un- 
certainty of our astrometric solution by adding in quadra- 
ture CTfit and the uncertainty a tI — ^/3/N s ctgsc = 0'.'12 
with which we can register our field on the GSC2 refer- 
ence frame (N s — 20 is the number of reference stars, and 
\/3 accounts for the free parameters in the astrometric fit; 
e.g., Lattanzi et al. 1997). The uncertainty on the reference 
star centroids is below O'.'Ol and was neglected. After ac- 
counting for the 0'.'15 mean systematic uncertainty on the 
tie of the GSC2 to the International Celestial Reference 
Frame (Lasker et al. 2008), the overall uncertainty of our 
absolute astrometry is O'.^S (0'.'19 per coordinate). 

The ATCA and Chandra positions of PSR J0108-1431 
are shown in Fig. 1, overlayed on the FORS1 U band 
image after astrometric recalibration. As seen from Fig. 
1, the backward proper motion extrapolation nicely fits 
the position of the object tentatively detected at the edge 
of the field galaxy. This is more quantitatively shown in 

1 http:/ /star-www. dur.ac.uk/j5draper/gaia/gaia.html 

2 http:/ /star-www.rl.ac.uk/Software/software.htm 
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Fig. 2. Proper motion measurement in right ascension (left) and declination (right) computed from the optical ( VLT), 
radio (ATCA), and X-rays (Chandra) coordinates of PSR J0108— 1431. The displacements are computed with respect 
to the first epoch VLT position of the counterpart candidate. 



as much as possible the signal-to-noise ratio, photometry 
was computed using customized aperture diameters of 0'.'6 
and 0'.'7, i.e. about the size of the image PSF in the B and 
U-band, respectively. Due to the uncertainty in the de- 
termination of the object centroiding (see §2.1), measure- 
ments were iterated centering the aperture at different po- 
sitions around the best centroiding estimate. Results were 
found consistent within the measurement uncertainties. 
Aperture correction was then applied using as a reference 
the fitted growing curves of well-suited reference stars se- 
lected in the field. Airmass correction was applied using 
the Paranal extinction coefficients measured with FORS1. 
Photometric calibration was applied using the mean, ex- 
tinction corrected, night zero point^f). 

We obtained U = 26.4 ± 0.3 and B w 27.9. As the 
higher background in the B band caused by the close-by 
galaxy makes our measurement very uncertain, we will use 
a fiducial error ±0.5 for the B magnitude. The object is 
undetected in the V band. We have remeasured the flux 
upper limit at the location of the candidate counterpart 
and obtained V > 27.8 at a 3a level, consistent with the 
limit reported by Mignani et al. (2003). The derived spec- 
tral fluxes are shown in Fig. 3. Since the Nh = 7.3 x 10 19 
cm~ 2 , estimated from the pulsar's dispersion measure, im- 
plies a reddening of w 0.05 magnitudes, i.e. much smaller 
than our photometric errors, we neglected the interstellar 
extinction correction. 

3. Discussion 

With only one fairly accurate measurement, our photom- 
etry does not allow us to determine the spectral slope for 
the candidate counterpart. The measured optical fluxes 



Fig. 2, where we compare the relative offsets between 
the coordinates of the optical source and those of PSR 
J0108-1431 measured by ATCA and Chandra. The coor- 
dinates of the optical source, computed with our astromet- 
ric solution, are a = 01 h 08 m 08?301 and 5 = -14°31'49"15, 
(J2000), with the la uncertainty of 0'.'27 in each coordi- 
nate. This value includes the overall uncertainty of our ab- 
solute astrometry and the object centroiding error, added 
in quadrature. Since it is practically impossible to fit a 
PSF to the object profile, we assumed a conservative es- 
timate of ±1 pixel (±0'/2) for the centroiding error. 

The yearly displacement computed using all the avail- 
able coordinates yields the proper motion of [i a — 100±41 
mas yr -1 , \ig — — 184±49 mas yr _1 in right ascension and 
declination, respectively. These proper motion values are 
in excellent agreement with those reported by Pavlov et 
al. (2008). Thus, based on the agreement between its posi- 
tion and the proper motion extrapolation, we propose the 
object originally identified in the VLT images of Mignani 
et al. (2003) as the candidate optical counterpart of PSR 
J0108-1431. 

2.2. Photometry 

We carefully remeasured the flux of the candidate coun- 
terpart. The object is fairly well detected in the U band 
(~ 5a), but it is only marginally detected in the B band 
(~ 3a). Due to the difficulties in fitting a model PSF to 
the object profile, we measured its flux through aperture 
photometry. Sky background was sampled in a number of 
regions selected within a radius of ~ 4" from the object 
(far enough from the field galaxy and from other relatively 
bright objects in the vicinity), with sizes 2 to 4 times larger 
than the FWHM of the image PSF. In order to increase 



3 http:/ /www.eso.org/observing/dfo/quality/FORSl/qc/qcl.html 
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are compatible with a thermal (Rayleigh-Jeans) spectrum 
emitted from the bulk of the neutron star surface with the 
(brightness) temperature T = (7-10) x 10 4 (d 130 /i? 13 ) 2 K, 
where i?i 3 is the apparent neutron star radius in units 
of 13 km. An example of the blackbody spectrum consis- 
tent with the optical fluxes of the candidate counterpart 
is shown in Fig. 3. 

The U and B band fluxes and the V band upper limit 
cannot be described by a flat or decreasing power-law 
spectrum, F v oc v~ a with a > 0, such as expected for 
magnetospheric emission. It could be described by an in- 
creasing power-law, but power-laws with a < have not 
been observed so far in the optical- UV spectra of rotation- 
powered neutron stars (see, e.g., Mignani et al. 2007). 
Thus, it is unlikely that the optical emission is purely 
magnetospheric. Of course, we cannot rule out that a mag- 
netospheric spectral component is present at a level below 
the V band limit (Fig. 3). 

The optical spectral fluxes are inconsistent with the 
extrapolations of the power-law and blackbody X-ray fits 
(Fig. 3). In the thermal interpretation, the X-ray emission 
comes from small, hot polar caps, whose contribution to 
the optical spectrum is negligible, but the optical emission 
may come from the bulk of the neutron star surface, which 
is too cold to be seen in X-rays. If the X-ray spectrum is 
predominantly magnetospheric, then the optical through 
X-ray spectrum of the magnetospheric emission cannot be 
described by a simple power-law, but it could be described 
by a model in which the spectral slopes are different in the 
X-rays and optical (e.g., ax = 1.2, ao = 0; see Fig. 3 for 
an example). 

If further observations confirm that the object is the 
pulsar counterpart with a predominantly thermal optical- 
UV spectrum (as suggested by Fig. 3), this would have 
important implications for understanding the thermal evo- 
lution of old neutron stars and the processes in their in- 
teriors. Current neutron star cooling models are rather 
uncertain for ages greater than ~1 Myr because the effect 
of possible (re)heating processes is unclear. Such processes 
include the Joule heating, heating caused by readjustment 
of the equilibrium neutron star structure in the course 
of its evolution, or heating caused by the effects of su- 
perfluidity in the neutron star interiors (see Fernandez & 
Reisenegger 2005, and Schaab et al. 1999 for a review). 
Evidence for reheating has been found in the much older 
recycled pulsar J0437— 4715, which showed a brightness 
temperature of ~ 1 x 10 5 K (Kargaltsev et al. 2004). The 
measurement of a similar temperature in PSR J0108— 1431 
would mean that some heating processes operate in slowly 
rotating, non-recycled pulsars. (We should note, however, 
that knowing the distance to the pulsar would be nec- 
essary to measure the surface temperature from the flux 
measurements in the Rayleigh-Jeans tail, so that measur- 
ing the parallax of PSR J0108-1431 is badly needed.) 

If further observations find a non-thermal component 
in the optical spectrum, this would be important for 
understanding the long-term evolution of the rotation- 
powered magnetospheric emission. Recent observations in- 
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Fig. 3. Results of the optical photometry of the PSR 
J0108— 1431 candidate counterpart (U and B band fluxes 
and V band upper limit) and the X-ray spectrum of PSR 
J0108-1431 measured with Chandra in the 0.3-8 keV 
band (Pavlov et al. 2008). The fits of the X-ray spectrum 
with the power-law and blackbody models are shown, with 
their extrapolations to lower energies. The hatched ar- 
eas correspond to la uncertainties of the X-ray fits. The 
dashed curve shows a blackbody spectrum with the tem- 
perature of 9 x 10 4 K (for R = 13 km, d = 130 pc). 
The dotted line shows an example of a flat power-law ad- 
ditional component in the optical-UV band, compatible 
with the measured V band upper limit. 

dicate that even very old pulsars are still efficient magne- 
tospheric emitters, with an even higher fraction of their 
spin-down power emitted in the X-ray and optical bands 
than in younger pulsars (Kargaltsev et al. 2006; Zharikov 
et al. 2006; Pavlov et al. 2008). In particular, in the op- 
tical/UV this is demonstrated by the cases of the pulsars 
B1929+10 (Mignani et al. 2002) and B0950+08 (Zharikov 
et al. 2004), which exhibit flat optical spectra with opti- 
cal efficiencies r\o = Lo/E <~ 1 x 10~ 6 and ~ 4x 10~ 6 , 
respectively (Zavlin & Pavlov 2004). If at least part of the 
detected emission from the pulsar counterpart candidate 
were of a magnetospheric origin, then the optical efficiency 
(e.g., rjo ~ 2 x 10~ 5 di 30 for the flat spectrum shown by 
the dotted line in Fig. 3) would be rather high compared 
to the younger pulsars B1929+10 and B0950+08. On the 
other hand, the (distance-independent) ratio of the X-ray 
and optical efficiencies (rjo/vx ~ 0.005 for the above ex- 
ample) is within the range found for other pulsars detected 
in the optical (see Fig. 6 of Zavlin & Pavlov 2004). 

4. Summary 

By re-analyzing the VLT images of Mignani et al. (2003) 
we have identified a posteriori a possible optical counter- 
part to PSR J0108— 1431, based on its positional consis- 
tency with the backward pulsar proper motion extrapo- 
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lation. If our identification is correct, this would be the 
tenth (or eleventh, if the detection of PSR Bl 133+16 by 
Zharikov et al. 2008 is confirmed) rotation-powered pulsar 
detected in the optical/UV. In any case, PSR J0108-1431 
would be the oldest non-recycled radio pulsar identified at 
optical wavelengths. As almost 8 years have passed since 
the reported observations, the optical identification could 
be easily confirmed by proper motion measurement of the 
candidate counterpart in deep optical/UV observations. 
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